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SUMMARY 

An optimal control solution is obtained for the descent and landing 
of a helicopter after the loss of power in level flight. The model considers 
the helicopter vertical velccity, horizontal velocity, and rotor speed; 
and it Includes representations of ground effect, rotor Inflow time leig, 
pilot reaction time, rotor stall, and the Induced velocity curve in the 
vortex ring state. The control (rotor thrust magnitude and direction) 
req,uired to minimize the vertical and horizontal velocity at contact with 
the ground is obtained using nonlinear optimal control theory. It is found 
that the optimal descent after power loss in hover is a purely vertical 
flight path. Good correlation, even quantitatively, is found between the 
calculations and (non-optimal) flight test results. The optimal control 
solution is thus a consistent and accurate method i >r comparing and evaluating 
the power-off descent characteristics of various he . <.copter designs . 


INTRODUCTION 

Good autorotation characteristics diirlng descent after power loss 
are essential for a useful and safe helicopter design. While it is known 
that the helicopter rotor has a minimum descent rate in vertical autorotation 
about the same as a parachute of equal size, there are other questions which 
require consideration. First, the helicopter rotor is a rather small parachute, 
so the ideal descent rate can be fairly high (V * 1.16 (T/a)^ m/sec, where 
T/a is the disk loading in kg/m ). This high basic descent rate increas'^s 
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the Importance of other parameters In the power-off landing maneuver. 

Secondly, it Is necessary to fly the helicopter in a manner to achieve 
the least descent rate, and most importantly to flare near the ground so 
that the helicopter touches down with vertical and horizontal velocities 
as nearly zero as possible. 

It is desirable to have in the preliminary design process a means 
of assessing the Influence of basic parameters on the helicopter autorotation 
characteristics. A number of elementary autorotation indices have been 
developed, generally based on the ratio of the rotor kinetic energy to 
the helicopter power required, KE/P (KE is the energy available during 
the descent, and P Is the rate of energy decrease just after the lose of 
engine power, thus a high ratio of KB/P is desired). The problem Is more 
complex really, with many parameters of the helicopter design influencing 
the autorotation characteristics. A difficulty lies in the necessity for 
flying the helicopter to the ground, which requires the choice of a control 
schedure. A poor choice for the helicopter control can easily result in 
an unacceptable landing, thus obscuring the influence of the design parameters. 
Therefore this report considers the use of nonlinear optimal control theory 
to establish the best control schedule, and thereby eliminate 
the influence of the control choice. The result is a consistent method 
for comparing and evaluating the power-off landing characteristics of 
various helicopter designs. 



EQUATIONS OF MOTION 


The optimization problem to be formulated is to find the control after 
powe.' loss to arrive at the ground with minimum veloclcy, given the 
helicopter Initial altitude h^, flight state, and basic parameters. The 
helicopter Is assumed to be in equilibrium level flight at the Instant of power 
loss, with rotational speed , rotor loading and forward 

speed . The basic parameters of the helicopter design include the Lock 
number ^ , the rotor radius R, and the solidity ratio . The aircraft 
position is defined by the coordinates h and t ., respectively vertical and 
horizontal (see figure l). It is convenient to measure h downward, so 
h * 0 at the initial altitude and h * h^ at the ground. 

The optimal control problem is best solved using an independent 
variable which has a fixed end point. Thus the Independent variable for 
the present problem must be tne haight h rather than time, since the 
arrival at the ground is defined by h « at an unknown time. The change 
of variables is accomplished using d( )/dt « h d( )/dh, or ( )* » A ( ) . 

The numerical integration which is required to solve the problem is still 
best done with respect to time however (see below). 


For the control variables, the magnitude and direction of the 
rotor thrust are used, specifically the thrust coefficient and the 
angle of the thrust vector to the vertical o< (see figure l)« It is 
convenient to express the problem in terms of the vertical and horizontal 
components of C,p, eose< and 6^ - C,p sino< respectively. The 

collective pitch control required to obtain this thrust nay be then obtained 
from the blade elemer.l theory emparession 


©TS ~ 
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It is not possible to obtain the longitudinal cyclic control from e< without 
considering the helicopter pitoh attitude and the rotor flapping alsoi the 
primary interest here is in the flight path anyway. 
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The equations of motion considered to describe the helicopter descent 
after power loss are those for vertical descent velocity h, horizontal 
velocity X, rotor speed and Induced velocity v. A separate differential 
eqtiatlon Is used for the Induced velocity partly to allow consideration of 
a time lag In the Inflow response, and partly to simplify the incorporation 
of the Inflow curve (including groiind effect) In the model. Vertical 
force equilibrium (see figure 1) gives i 


Mh “ W - T cosav + D slnd 


or since W > Mg and h 


• v 

Jl. * 


V ' W vsj ^ 




Horizontal force equilibrium gives 


Mx “ T slnw - D cosd 


I ( 




The helicopter parasite drag will be defined by an equivalent area f , 
such that D “ -I V^fj then D sin© = -|-j^ h(h^+ x^)^f, and D cosS =* 
■|■;^x(h^ + x^)^f (see figure l). Rotor torque equilibrium after the loss 


of engine power Isi 


h-^ - - ft 


XI — ~ 




where Q Is the rotor aerodynamic torque, given by 

(reference l). Here is the rotor stall limit (with n a large number, 

IS 8 

®*8* ■ 20, so the profile torque greatly increases when the loading is above 

D 

C^er^), The rcior advance ratio ^ and inflow ratio are given by 









X e»So< — f~ 
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The rotor Induced velocity v is given hy a differential equation which 
Includes a time lag t 

TV + V = Kvj^fjf^ 

or 

"*• Vfr 

'c 

The steady state solution is thus v « vcv^fjfj,. 

Here v^ - T/2^A, »c is an empirical factor (typically ^ “ l*15)t 
f. Is the Inflow curve, and f- is the effect of the ground. 

For the inflow curve, the following expression is usedj 
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where the parameters x and y are defined by 
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The first expression for fj Is the usual moraentua theory result (ref. l){ the 
second expression is an empirical approximation for the vortex ring 
state (where the momentum theory hreaks down). The region of rou^ness 
in the vortex ring state Is defined approximately by (2x+2) +y <1. 

To account for ground effectt the following expression is used (ref. 2): 


Here z is the rotor height above the ground, z-h - h + z (z is the 

o o o 

rotor height for the aircraft on the ground); and £ is the angle of the 
wake to the ground. 
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DIMENSIONLESS EQUATIONS 


Now the equations will be made dimensionless using the quantities 
• and R • The four degrees of freedom are defined as follows t 

d » h/Xl,R 
e - x/HoR 

jt - v/SlJi 

end the control variables are cos&< and «• sln«K , The 

four dlffesrentlal equations are them 

^ 3: ( I - ^ ') 

® A Cto ^ 

3: ( ^ c y 

^ ^ ^ Cx. 
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Frr and f^ In the Inflow equation, the following quantities are required i 

eCa — 
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Finally, the initial conditione (for level flight) are d •■ 0, e - * 


<^ ■ 1, and ■ J?o ■ v^/Si»R at h * 0. 


CRIISRION 


The optiaization problem ia to arrive at the ground with minimum 
vertical and horizontal velocity. Thus a quadratic cost function of 
the following form is usedi 

2 . „ »Z 


J - i ( hf + «x ^f ^ 
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where and are the velocities at the ground (h • h^), and Is the 
weighting function of horizontal velocity relative to vertical velocity. 
Sow since h • ^ (ih^) ■ g(l - Tcos«/h + DslnB/w), there follows 
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So an equivalent cost function 1st 




In terms of the dimensionless quantities then. 
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The control prohlen Is to find C and C as a function of h to minimize J, 

8 X 

subject to t«M constraints defined by the differential equations above. 
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HONIJMBAR OPTIMAL OOMUfOL 

Consider a system defined by the nonlinear diffe7.*ential equation 
X • a(x, u, h), where x is the state vector i u is the control vector, 
and h is the independent variable; and a cost function J = b(x, ui h) dh. 
It is assumed that 'Uie initial conditions x(h^) are given, and that h^ 
and hj are fixed. The optimal control problem 1^ to find u(h) to minimize J. 
The solution (see reference 3) is defined by the following set of equations i 


X = ov 


•ay 

f - 


- Y 


hW / 6 <K ^ 

Sc ~~ hC ^ 

with boundary conditions x(hj^) • and p(h^.) ■ 0, 

For simple problems, -toe equation ■ 0 is solved directly 

for u as a function of 'x, p, and h; and u is substituted into the first 
two equations. The differential equations are then integrated, using the 
boundary conditions to eliminate integration constants. Then the solution 
for is the optimal path, and p gives the optimal control law u(h). 


The present problem is too complex for such a procedure, so a 
steepest descent algorithm is used to solve the two point boundary value 
pi<obl«B (reference 3). A cycle in the algorithm consists of the following 
steps. A current estimate of the optimal control law, u , is available. 
The differential equation x^ - a is integrated forward from hj^ to h^ using 
u'**' and the Initial consitlons on x. Next the differential equation 
p _ h is integrated backward from h^ to using x, u*^*" , and 

the final conditions on p. Finally, ( ^ is evaluated using x, p*, 

and u^'*' I and the control is incremented by 


tA — *A 


Si 


where ^ is a step sise, chosen (by trial and error) such that the solution 






converges fast enough without overshooting* This process is repeated 
until the solution converges to the optimum, as indicated by the cost 
J approaching a minimxim* Such a steepest descent procedure has the advantage 
of not being sensitive to the initial guess for the control} the convergence 
slows down as the minlmim is appu.^achei however. 


OPTIMAL CONTROL PROBLEM 

The optimal control problem for the power“Off descent and landing 
of a helicopter is obtained by applying nonlinear optimal control theory 
to the dimensionless equations given above. The formulation using h as the 
independent variable is required since the final height is specified , rather 
than the final time. It is still preferable to do the actual numerical 
integration using time as the independent variable however, to exlminate 
the singularity which occurs at d ■ 0 (such as at the start of the maneuver) 
is h is 'O’.e inov'pendent variable. Therefore after the differential equation 
for is opbtained, the coordinate txansfomatlon is maide back to t, using 

d( )'*^ “ ( y • is also necessary then to integrate h ** iCiji to obtain 
the proper variable h (and also x • ei\R to obtain x( t) ) . The resulting system 
of equations for the optlMl control problem is then ss follows. 

y - a j J ) 
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and the natrlcea B and C are defined as follows. 
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These equations are solved by the steepest descent method outlined 
above. The pilot reaction time tp is included by constraining the collective 
pitch to have the initial value for t < tp. Then the control is given by 

Cl T- — - 1 ■ 

2 

Assuming o< is unchanged for t < tp, then /Czo “ (if/A/^r)/CTo. 

Typically tp ^ .75 sec (althouf^ handling qualities fipeclficat^ons may 

require the use of a larger value). 


IfiSCBIIT FROM HOVER 

Consider the caise of optimal descent after power loss in hover, 
hence with initial condition 0. The solution of the above equations 
will be shown to be ■ 0 and e * 0 (o< =0 and x » O). 

Assuming ■ e ■ 0, it follows that 0, y ■* 0, ^fj/^y * 0 
(l.e. /V- 0 at “ O), and f^ is a function of h only. The 

differential equation for e becomes e * 0, with solution e * constant ** 0 
(xislng the initial condition e “ >Uo " O). The differential equation above for 
p^ becomes then pj ■ /^o^^e* ^ solution 

Pe ” Pe^^^ t«o U 

or p s 0 since the final conditions on give p (h ) = 0. Then bH/<^C = 
g (o^ydCf )p 3 0, as required for the optimal solution. The remaining problem 
has then only three degrees of freedom (d, oj , and ), one control variable 
* C^), and threee components of p (p^, p^ , and Pj^ ). While eliminating 
e, C , and p from the problem is a significant simplification, it is still 
necessary to integrate numerically and i cerate to find the optimal solution. 

Thus the optimal contzttl solution for descent from howver after power 
loss is a purely vertical flight path. The same conclusion was reached from 
the flight tests reqorted in reference 4, although in practice a small amount 
of forward speed is required, both to avoid the vortex ring state during 
flare and to keep the landing point in sight. 
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RESULTS AND DISCUSSION 


The optimal descent of a helicopter after power loss has been 
calculated for a nvanber of cases. Because It is found both by flight 
tests and from calculations that an initial foiward velocity greatly 
Improves the autorotation characteristics, res\dts are given here only 
for descent from power loss in hover (which as found above involves 
a purely vertical flight path). The helicopter considered is one for 
which flight test data are available. Three values of the rotor Lock 
number are considered, from ■ 4.5 (the standard rotor) to ^ * 2.6 
(a rotor with heavier blades, investigated specifically for better 
autorotation characteristics). Figure 2 shows xhe vertical velocity at 
the instant of contact with the ground, after optimal descent from power 
loss in hover at varous altitudes. These calculations are in agreement 
with the flight test results. Specifically, it was found in reference 4 
that the autorotation characteristics greatly improved for the heavier 
rotor (the autorotation characteristics for the helicopter with )( * 4.5 
are poor in this altitude range, while t}» characteristics of the helicopter 
with ^ « 2.6 were found to be very good)i and the critical helj^t with 
■ 2.6 was about h^ ■ 30m. 

Figure 3 presents in detail the optimal solution for power-off 
descent from hover at an altitude of h^ *■ 30a. Figure 3(a) gives the 
collective pitch control required as a function of tine (note that a pilot 
reaction time of .75 sec is used). These res\ilts again agree generally 
with reference 4 , which found in flight tests that the collective should be 
dropped immediately, followed by a gradual increase for flare (in figure 
3(a) the flare begins when the helicopter is about 20m above the ground). 
Figure 3(b) shows the rotor as a function of altitude (the control 

variable actually used In the solution procedure). A stall limit of ^ 

^ ” .15 was used, which results In a leveling off of the control 
Just before It reaches that value t the hlc^ torque due to rotor stall 
greatly slows down the rotor, and thus values of above stall are 

not called for until Just before ground contact. Figure 3(e) shows the 
helicopter vertical load factor (n„ ■ t !3/^vo)i 

Si 


- 15 - 


figure 3(d) shoKs the rotor speed, as a fraction of the Initial value i and 
figure 3(e) shows the vertical descent velocity (the rotor Is operating 
In the windmill brake state when the velocity Is Increasing:, and In the 
vortex ring state at the end of the maneuver when the velocity is decreasing) 
Finally, figure 3(f) presents the flight path for the optimal descent; 
note that the principal influence of Lock niaiher Is on the final portion of 
the flare for these caises, where the extra kinetic energy in the heavier 
rotor allows a greater reduction In velocity. 

Figures 4 and 5 *present a comparison between flight test results 
(unpublished data from the program reported In reference 4), and the 
calculated optimal power-off descent from hover. An optimal flight path 
was not flown in the tests of course, and In addition some forward speed 
and cyclic flare were involved. Sven so, the correlation Is qualitatively 
good. The most important discrepancy is that the rotor speed in the flight 
tests loss not Initially decrease as fast as in the calculations. Sxaminlng 
the flight test data, it is found however that the engine torque does not 
decrease to zero immediately after the throttle chop; in fact, vhe torque 
remains above of full power until the helicopter has descei^ed about 
15m. Adding to the analytical model an exponential lag in the engine 
power drop greatly improves the correlation in figures 4 and 5* While 
again one should not look for too much correlation here, the lower vertical 
load factor in the flight tests during ^e collective flare suggests that 
the ground effect might be stronger than was used in the model, perhaps 
due to the helicopter vertical velocity (the cyclic flare may be Influencing 
the meastired n^ also however). 
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OONCLUDIIIG REMAHES 


An optlaal control solution has hesn obtained for the descent 
and landing of a helicopter after power loss. A comparison with flight 
test results shows sufficient correlation, even quantitatively, to verify 
the basic features of the model* The Influences of parameters such as 
altitude and Lock number are correctly given, and the proper characteristics 
of the control technique are obtained. This model should thus prove to be 
a useful tool for evaluating and comparing the power-off landing 
characteristics of various helicopter designs. 
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Figure 1 Dafinltlon of helicopter position (x and h), velocity 
(V and d )• and forces (V, D, T, and ). 
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Figure 3 Optlnal power-off descent from hover at altitude 

= 30 m, for three values o"f* rotor Lock number ^ 
(R “ 5*38m, nR» 199 m/sec, “ .063, and 

.046) 
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Compiete power loss 

With engine power lag 

O Flight test 





Comparison between flight test data and optimal power-off descent 
from hover at altitude h = 38m. (R - 5 . 38 m, r2.R = 199 m/sec, 
V'^o " -066,^- .048,° t - 2.6) 
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